Witnessing non-classicality of a quantum oscillator state by coupling it to a qubit 
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We propose a new witness operation for the non-classical character of a harmonic oscillator state. 
The method does not require state reconstruction. For all harmonic oscillator states that are clas- 
sical, a bound is established for the evolution of a qubit which is coupled to the oscillator. Any 
violation of the bound can be rigorously attributed to the non-classical character of the initial 
oscillator state. 
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I. INTRODUCTION 

Quantum mechanics is fundamentally different from 
classical mechanics. Interference of probability ampli- 
tudes, superposition of states, uncertainty relations be- 
tween canonically conjugate variables, etc., are essential 
quantum phenomena that are not present within the clas- 
sical theory. If a state of a system exhibits any such 
intrinsically quantum feature, the state is called non- 
classical [1]. 

For a given physical system, quantitatively categoriz- 
ing the states into classical and non-classical is usually 
challenging. In this report, we are interested in catego- 
rizing the states of a harmonic oscillator. Any harmonic 
oscillator state can be written in the coherent state diag- 
onal representation 0, H| : 



(1) 



p = / d z aP(a)\a)(a\, 



where |a) is a coherent state. Within the fundamen- 
tal limits imposed by the uncertainty relation between 
the position and the momentum, a coherent state cor- 
responds as closely as possible to a classical harmonic 
oscillator of a definite complex amplitude. For this rea- 
son, a coherent state can be considered to be classical. A 
natural definition of classicality based on this observation 
was introduced by Glauber Q. He proposed that if P(a) 
is a valid probability measure, the state p can be thought 
of as a statistical mixture of various classical states and 
thus is classical itself. On the other hand, if P(a) is not 
a valid probability measure, the state p is non-classical. 
We will adopt this definition of non-classicality of a har- 
monic oscillator state. 

There are many ways of checking whether the P dis- 
tribution corresponding to a given oscillator state fails 
to be a valid probability measure. For example, one can 
perform a complete state tomography to find out what 
P is [H HI ■ Knowing the state, one can then derive the 
P distribution. Such an approach, although possible in 
principle, is difficult to carry out in practice because the 



inverse relation connecting the distribution to the state 
involves an integral with an exponentially diverging fac- 
tor m: 



P(a) 



d 2 7e' 7 ' (— ~f\p\~f)e~ 



(2) 



Because of the diverging term, any experimental error 
gets exponentially enhanced Q ■ For this reason, instead 
of performing a full state tomography, one generally looks 
for non-classicality witnesses that signify P as an invalid 
probability measure. Some of these witnesses are anti- 
bunching Q, sub-Poissonian statistics [H, [lfj, quadra- 
ture squeezing [llj], slower decay of the characteristic 
function of the rotated quadrature distribution than the 
characteristic function of the ground state [13] , negative 
Wigner distribution of the oscillator state Q, violation 
of Bochner's criterion for the existence of a valid positive 
semidefinite characteristic function of the P distribution 

[nm, etc - 



II. WITNESSING NON-CLASSICALITY 
THROUGH A QUBIT COUPLED TO THE 
OSCILLATOR 



In this report, we propose another observable signa- 
ture of non-classicality of an oscillator state. We longi- 
tudinally couple the oscillator to a qubit with the joint 
Hamiltonian given by: 



H = n — er 2 
2 



hue 
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(3) 
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where ujq and to are the qubit and oscillator frequencies 
and (3 is a dimensionless, constant coupling parameter. 
The oscillator operators, a and a\ are the annihilation 
and creation operators respectively and the qubit oper- 
ator, er 2 , is a Pauli matrix. This Hamiltonian is differ- 
ent from the Rabi Hamiltonian where the qubit-oscillator 
coupling is through a x . 

Because of the qubit-oscillator interaction, it is possi- 
ble to learn about the initial state of the oscillator by fol- 
lowing the evolution of the qubit. The method of recon- 
structing an arbitrary oscillator state by looking into the 
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dynamics of an interacting few level system was proposed 
in Later, using an interacting qubit, the method 

of state reconstruction was developed in areas of cavity 
QED [H-SIl, circuit QED [H, M, trapped ions [HHl] 
and nano-mechanical resonators |27l. l28j . 

It was shown by Tufarelli, et al., that by using a Hamil- 
tonian that is similar to Eq. © , one can reconstruct the 
entire state of the oscillator 28] . This state reconstruc- 
tion method is experimentally extremely challenging as 
it requires a predefined, deterministic modulation of the 
coupling strength in time. Also in the reconstruction 
procedure, it is necessary to be able to change the am- 
plitude of modulation in different experimental runs. In 
the present report, we propose a different method of wit- 
nessing the oscillator state non-classicality through the 
qubit dynamics that can be carried out using a constant 
coupling parameter, j3. 

We assume that the qubit and the oscillator states are 
initially separable and their joint state is given by: 

Pqo = Pq® Po: 

y?J) 1 ! ( ° ( ) 0) )«/d 2 aP(a)|a)(a|. (4) 

The oscillator state, p , is the state whose non-classicality 
we want to investigate. The rows and columns of the 
qubit state are defined in the basis: |+) and — ), where 
<r z |±) = ±|±). The parameters, z and w, defining the 
matrix elements of the qubit are related to the expec- 
tation values of the Pauli matrices: z = ((<r z ) + l)/2, 

w = ((o-*) -*K))/ 2 - 

Using the Hamiltonian, Eq. (|3|), the time evolution 
of the joint state, p qo , can be found. The qubit state 
can be evaluated by tracing over the oscillator degrees 
of freedom from the time evolved joint qubit-oscillator 
state: 



Pq {t) = Tr {e- lHt ' h Pqo e lHt / h }. 



(5) 



The diagonal terms of the qubit density matrix, writ- 
ten in the a z eigenbasis, do not change in time, i.e. 
z(t) = z(0). This is because a z is a constant of motion: 
[a z ,H] = 0. For the off-diagonal terms, we find: 



>(t) 



—ILJQt 



g(t)w(0), 



(6) 



where g(t) depends only on the oscillator degrees of free- 
dom and is defined as 

g(t) =Tr {e- lff +*/ ri p e 4H - t/?1 }, 
d 2 aP(a) 

xTY {e' lH + t/h \a){a\e lH - t/h } 1 (7) 



where 



H± = ^(a t a±/3(a + a t )) 



(8) 



The H± operators correspond to the Hamiltonians of dis- 
placed harmonic oscillators. So, one can evaluate g{t) 
analytically to get: 



x e 



-4i l 3(ae~ i "" /2 +a*e i '"* /2 ) sin ^ 

where we have defined 

f(a,t) = e -^(«e- tot/2 +^ t/2 )^¥ 

and 

W(t) = ( d 2 aP(a)f(a,t). 



(9) 
(10) 

(11) 
(12) 



If we look at the absolute value of the function g(t), 
we get from Eq. ([9]): 



\g(t)\ 



< e 



-8£ 2 sin 2 (^) 
-8/3 2 sin 2 (if ) 



d 2 a P(a) f (a, t) 



P(a)f(a,t) 



(13) 



If P{a) is a valid probability measure, we can write 



P(a)f(a,t) =P(a) f(a,t) 



(14) 



Putting Eq. in Eq. (TT3]) and using the fact that 

\f(a,t)\ = 1 and J d 2 aP(a) = 1, we get: 



\g(t)\ < e - 8 ^ si " 2 (^) j d 2 aP(a), 

= g-S^sin 2 ^ ) 



(15) 



Using this upper bound for g(t) in the expression for the 
off diagonal element of the qubit density matrix, Eq. ^ , 
we get: 



|IU(t)| = e 8 ^ sin2 (^) 



w(t) 



w(0) 



< 1. 



(16) 



Note that w(t) and thus IW^)! can be experimentally 
measured by measuring the expectation value of the 
qubit observables, a x and (7 y , CIS Ci function of time. 

Inequality (|16j) is the main result of the paper and 
we call it the non-classicality witness inequality (NCWI). 
The inequality states that if the oscillator state is such 
that the associated P distribution is a proper probability 
measure implying Eq. (|14p to be correct, |VF(i)| is always 
bounded from above by unity. Thus if |lU(i)| > 1, we 
know for sure that the initial oscillator state, po, is non- 
classical. 

Two remarks are in order at this point. First, it can be 
shown that W(t) is related to the characteristic function 
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of the P distribution [29f. This suggests that Eq. ([16]) 
is related to the criterion for observing non-classicality 
through measuring the probability distribution of rotated 
quadratures [l2[ ■ The second remark is just that a viola- 
tion of the NCWI is only a sufficient but not a necessary 
condition for determining whether a state is non-classical, 
i.e., a state can be non-classical and still not violate 
the witness inequality (see the examples of Fock states, 
vacuum subtracted thermal states and Schrodinger Cat 
states in the next section). 



III. EXAMPLES 

To illustrate the above statements, we calculate | W(t) | 
for various initial oscillator states. 

Coherent state: The P distribution for a coherent 
state, say |ao), is P co h{a) = 5 2 (a— ao). Corresponding to 
Pcoh( a )i we have |W co h(£)| = 1- Thus, no coherent state 
ever violates the NCWI. This is required for consistency 
because by our definition a coherent state is considered 
to be a classical state. 

Thermal state: For a thermal state with mean excita- 
tion number n, the P distribution is given by: 



Pth(a) 



1 



■ exp (— \a\ 2 /n). 



(17) 



We see that Pth(ct) is a Gaussian probability distribu- 
tion and thus a thermal state with arbitrary n can be 
considered classical. Using P t h(a), we get 



\W th (t)\=e' lm02sin2 ^). 



(18) 



In agreement with our intuition that a thermal state is 
a classical state, we see that |Wt/,(i)| is always less than 
one and thus never violates the NCWI. 



Fock state: The P distribution for a Fock state, \N), 



is 



Pn{o) 



exp(aa*) / d 



32JV 



da N d 



*N 



5\a) 



(19) 



Because Pn (a) for N > involves derivatives of the delta 
function and is more singular than the delta function 
itself, each Fock state, except the ground state, is non- 
classical. Calculating Wjsr(t) using Ppf(a), we get: 



\W N (t)\ 



L N 16^ sin" 



uit 



(20) 



where Ljv(x) is a Laguerre polynomial. 

In Fig. [TJ we plot |Wjv(^)I f° r various values of N and 
for a given constant coupling strength, (3 — 0.5. Noting 
that |Wjv(t)| is a periodic function, we plot |Wjv(t)| only 
for a single period. We see from the figure that for some 
Fock states, e.g. for N = 1 and 10, the NCWI is violated. 
This agrees with the non-classical nature of these Fock 
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FIG. 1. (Color online) The function \W N (t)\ for = 0.5 and 
various values of N. We clearly see that |Wjv(t)| violates the 
NCWI indicating that Fock states are non-classical states. 



states. On the other hand, it is important to note that 
for any given coupling strength, not all Fock states will 
violate the NCWI. For example, in Fig. [TJ no violation is 
shown by |Wo(t)| and | W15 (i) | ■ However, this result does 
not suggest that some Fock states are classical. This 
point is further examined in the following examples of 
vacuum subtracted thermal state and superposition of 
classical states. By looking at the as ym ptotic expression 
of a Laguerre polynomial for large N [30( , it can be shown 
that for a given coupling strength, f3, the NCWI will not 
be violated for highly excited Fock states. 

Vacuum subtracted thermal state: We now consider the 
state 



Pth 



J22- N \N){N\. 



(21) 



N=l 



The above state corresponds to a thermal state with 
unit mean excitation number and from which the vac- 
uum state has been discarded. The P distribution for 
this state is 



PM = - C x P (-\a\ 2 )-S 2 (a). 



(22) 



We see that P th (a) has a negative measure at a 

a 





This implies that p' th is a non-classical state. Corre- 
sponding to P th {a), we have 



\Wl h (t)\ = |2e- 16 ^ sin2 (^) -1|. 



(23) 



Although P th (a) corresponds to a non-classical state, we 
see that |W t ' ft (i)| is always less than one and never vi- 
olates the witness inequality. This illustrates the fact 
that violation of the NCWI is not a necessary but only a 
sufficient condition for determining non-classicality [3l|. 

Superposition of classical states: Even though a co- 
herent state is considered to be classical, a state consist- 
ing of a superposition of coherent states might exhibit 
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FIG. 2. (Color online) The function \W ao (t)\ for /3 = 0.5 and 
various values of ao- For some values of ao, |W / sc (i)| violates 
the NCWI. This indicates the non-classical nature of these 
Schrodinger cat states. 



non-classicality. To understand this, let us examine the 
Schrodinger Cat state which consists of an equal super- 
position of two coherent states, one being the negative of 
the other: 



|* sc (ao)> 



K) + 1 - a ) 



2(1 



-2a 



°) 



(24) 



For simplicity, we have taken ao to be real. The P dis- 
tribution corresponding to \^ sc (ao)) is jT|: 

P sc (a) = N 2 \d 2 (a - a ) + 5 2 (a + a ) + e (H 2 H"ol 2 ) 



(« 



- aod «-e aod As 2 (a)\, (25) 



where l/Af 2 = 2(1 



-2a 



°). We see that for ao 7^ 0, 



P sc (a) contains infinitely high order derivatives of the 
delta function. This clearly suggests that a Schrodinger 
Cat state is a non-classical state if ao ^ 0. 

For the Schrodinger Cat state, W sc (t) can be calculated 
analytically to get: 



ao = 1 and 2. This confirms the non-classical nature of 
these states. On the other hand, for any given value of 
/?, we notice that not all Schrodinger Cat states violate 
the NCWI. This is evident from the plots corresponding 
to ao = and 5. 



IV. CONCLUDING REMARKS 

In this report, we use Glauber's criterion for non- 
classicality, which is the failure of an oscillator's coherent 
state P distribution to behave as a valid classical proba- 
bility distribution. We demonstrate a witness operation 
for non-classicality of this type. We show that one can 
imprint a signature of non-classicality onto the dynamics 
of an interacting qubit. The result is that all classical 
initial states of the oscillator lead to a fixed qubit bound, 
given in equation (ITB1) . In other words, the qubit oper- 
ator's time evolution, expressed via |W(f)|, must remain 
bounded from above if the oscillator's initial state was 
classical. A violation of the bound is a direct indication 
of the non-classical nature of the initial oscillator state. 
The method does not require state reconstruction, and 
monitoring of the qubit can be confined to a single period 
of the oscillator. A number of examples are presented to 
illustrate the behavior of the bound, and confirm intuitive 
expectation in special cases. This strategy for witnessing 
non-classicality is well suited for physical systems where 
strong longitudinal coupling can be achieved between a 
qubit and a single-mode of an oscillator (3^, IH[ • Because 
of the unavoidable interaction of the qubit-oscillator sys- 
tem with its environment in any experiment, it is im- 
portant to take into account the effect of noises on the 
evolution of |W(i)|. Under the Markovian approximation 
of system-environment interaction, one expects |W(£)| to 
decrease exponentially in time and this decay can be com- 
pensated for if it is experimentally possible to obtain the 
exponential decay constants [28j . 



|W ac (i)| = N 2 cos(4/3a sinwf) 

+ e- 2o? cosh f 8/3a sin 2 ^ 



(26) 
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